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Abstract 

We present a program to evaluate tree-level multi-gluon amplitudes with up to 
two of them off-shell. Furthermore, it evaluates squared amplitudes summed over 
colors and helicities for up to six external gluons. It employs both analytic expres¬ 
sions, obtained via BCFW recursion, and numerical BCFW recursion. It has been 
validated numerically with the help of an independent program employing numerical 
Dyson-Schwinger recursion. 


1 Introduction 

Factorization prescriptions are powerful tools to tame the eomplex calculations involving quan¬ 
tum chromodynamics (QCD) for scattering processes at collider experiments like at the Large 
Hadron Collider. They factorize eontributions to cross sections according to the seales involved, 
and/or aceording to universality and aeeessibility via perturbation theory. Many faetorization 
preseriptions are heuristic in nature, and some are proven, whieh means that perturbation theory 
and the treatment of possible singularities, along with the oeeurrence of large logarithms of ratios 
of seales, ean be dealt with in a systematie manner. 

Partonie seattering amplitudes form an essential ingredient in faetorized ealeulations of eross 
seetions for hadron collisions. Even though the partonie states are not physical, on-shell par¬ 
tonie amplitudes are well defined gauge invariant objeets of the gauge theory QCD, using its 
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Lagrangian and the Lehmann Symanzik Zimmermann reduetion formula. Factorization embeds 
the non-physical scattering amplitudes into physical cross sections. Recently, it has been shown 
that scattering amplitudes involving any number of ojf-shell external gluons can also be defined 
in a rigorous manner m. Such amplitudes are relevant in factorization prescriptions requiring 
off-shell initial-state partons, like High Energy Factorization (HEF) |l2l[3l. Recent developments 
and calculations involving such factorization prescriptions can be found in ll4]-[T^ 

Calculations employing collinear factorization, for which the scattering amplitudes are com¬ 
pletely on-shell, have been automated to the end for arbitrary processes, with essentially arbitrary 
multiplicities, and within essentially arbitrary models of quantum field theory lfT3l - [23 . By now, 
developments are heading at reaching this status to next-to-leading order in perturbation theory. 
This includes one-loop amplitudes and real-radiation contributions with all the complications 
arising due to mass singularities and the highly non-trivial phase space integration Il23l - l29]| . 

HEF requires partonic scattering amplitudes with off-shell initial state partons, and automa¬ 
tion of the calculation of these has not been established. Systematic formulations of their calcula¬ 
tion have been established ll30] - [34ll . In this paper, we present a program to numerically evaluate 
tree-level multi-gluon scattering amplitudes with up to two off-shell gluons as function of the 
gluon momenta, squared and summed over colors of all gluons, and summed over the spins of 
the on-shell gluons. It evaluates them via color-ordered helicity amplitudes that are calculated 
using the generalization of Britto-Cachazo-Feng-Witten (BCFW) recursion [[35l[36l|, described 
in [[T7]| . to include off-shell gluons. The program uses both hard-coded expression obtained via 
analytical BCFW recursion, and numerical BCFW recursion. Using the latter, color-ordered am¬ 
plitudes may be calculated to essentially arbitrary multiplicity. Squared amplitudes summed over 
colors and helicities are provided for up to six external gluons. 

This paper continues as follows: in Section [2] the amplitudes that the program calculates are 
defined. Section [3] explains how color is treated. Section |4] describes the usage of the program, 
and Section |5] introduces the program with which it was validated. Section finally, contains 
the summary. 

2 Definitions 

We consider the generally factorized formula for the gluonic contribution to a cross section in 
hadron collisions 


&(g*(ki)gnk2)^X) 


(r(hi(pi)H2(p2] ^ X) = d% Fi(ki) d4k2F2(l<2) 


4v/(krk2p-kfk2 


This formula is very general, and g* does not necessarily refer to an off-shell gluon. In collinear 
factorization, for example, we would have 
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where fi is the eollinear pdf for a hadron of type i. We include the factors establishing averaging 
over spins and colors in here. In the hybrid HEF If^ . for example, F 2 would be as above, 
while Fi would be given by 


Fi(ki) = 




d^kj 

In 


dxi 

0 


3"i(xi,kT, p.] 6^(ki -xi Pi -k- 


tJ , 


(3) 


where S^i (xi, kj, p) is the unintegrated gluon density. 

The symbol X stands for a partonic final state, for example a number of on-shell gluons. The 
partonic cross section & is given by 


er(g*(ki)g*(k2)^X) 


dcI)(ki,k 2 ^X) |,A(g*g*^X)|'0(X). 


(4) 


The phase space integration includes the summation over all color and spin degrees of freedom. 
The observable 0 turns the partonic final state into a physical final state, for example through 
a jet algorithm. This function also contains the necessary symmetry factors related to the final 
state. We concentrate on the amplitude —> X) from now on, for the case that X stands 

for a number of on-shell gluons. 

We will adopt the convention that momenta denoted by the letter p are always light-like, while 
momenta denoted by the letter k are not necessarily light-like. For the multi-gluon amplitudes 
we consider, this means that the initial-state momenta are denoted k'j^, k^, while the final-state 
momenta are denoted pj,..., pj{^. Momentum conservation is imposed as 


ki" + k2 + P3 H h pjj — 0 , 


(5) 


so the initial-state momenta have negative energy. 

Spin amplitudes with on-shell gluons depend on the momenta p{^ and polarization vectors ef 
associated with those gluons. On-shellness implies that for each on-shell gluon we have 


Pi'Pi = 0 and pi-£i = 0. 


( 6 ) 


Gauge invariance assures the Ward identity that any momentum proportional to p|^ may be added 
to ef without changing the amplitude: 


pf,£f;...) pf, + zp[^;...) Vz G C 


(7) 


Consequently, the amplitude only depends on the transverse part of ef, and one may consider 
helicity amplitudes instead, which depend on p[^ and the helicity At which takes the two possible 
values -F or —. 

Regarding off-shell gluons, the amplitude depends on their momenta kf, and on their “polar¬ 
ization vector” or direction pf, satisfying 


PfPt = 0 and pfkt = 0. 


( 8 ) 
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Within HEF, the direction is given by the momentum of one of the scattering hadrons, and the 
off-shell momentum is typically defined in terms of this and a transverse momentum via 




(9) 


Realize, however that given any momentum kf, one may construct an associated direction pf 
satisfying Eq. ([8]) (see Appendix lAl). For the amplitude, the notion of transverse momentum is 
arbitrary to a certain degree, since one can shift a fraction of pf to kj^: 


i<“ = <p“ + q" with 


= Xi — X 


Ki = K,i + ^Vi 


( 10 ) 


Interpreting this as a change of pf rather than xi, the amplitude scales homogeneously with this 
change: 


A( ...;k^^p^,. 

Xi 


X 


= ^yi(...;kr,pr; 

Xi 


( 11 ) 


3 Color treatment 


One issue regarding the color treatment as presented in ll^ has to be settled. There, the am¬ 
plitude with off-shell gluons is obtained by considering each of them as an auxiliary eikonal 
quark-anti-quark pair, carrying fundamental color indices. This situation is different from n], 
where each off-shell gluon carries a single adjoint color index, and the well-known color de¬ 
compositions, the one presented in [|^ in particular, hold manifestly. It is not a priori clear that 
these also hold in the formulation of If^ . In particular, one might expect that the so-called U(1 )- 
gluons would contribute in the color-flow representation, connecting two eikonal quark lines. In 
this section, we will argue that this is not the case, and that the representations of |[I1 and ll^ 
are indeed equivalent. 

This is essentially guaranteed by the (proven) observation, stated in eq.(40) and eq.(41) 
of [[33l, that the induced vertices of figure 4 in If3^ are traceless. This also guarantees the 
equivalence of these vertices with those defined in ll^ . Consider an amplitude in the color 
representation in which each gluon, be it on-shell or off-shell, is represented by a pair of fun¬ 
damental color indices. For the on-shell gluons this is achieved by contracting the amplitude 
with T“{^, where ai is the adjoint color index of the external on-shell gluon. Tracelessness of the 
induced vertices implies tracelessness with respect to the color indices of the auxiliary eikonal 
quark-anti-quark pair, so for each pair of fundamental color indices referring to a gluon in the 
amplitude, on-shell and off-shell, we have 

( 12 ) 


Each pair i^, ]i may refer to a gluon, on-shell or off-shell, or an ordinary quark-anti-quark pair. 
The relation above only holds if g refers to a gluon. The general formula for the color-flow 
decomposition of the amplitude is given by II40[I41II 
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where the sum is over all n! permutations a of (1,2,..., n]. Thanks to Eq. (fT^ . all partial 
amplitudes vanish if I = cr(l) for any I refering to a gluon. As a result, for the ease there are 
only gluons, the decomposition reduces to Il42ll 


= 

n "On 


L 




) ct (2) ] ct (3) 


]n 


)a(l) 


O-eSn-l 


(14) 


where the sum is now over only (n— 1)! permutations. The change in notation from to A(j just 
indicates that the labelling is not identical, e.g. A 1234 = 4 I 2341 etc.. We may return to the adjoint 
representation by contracting every fundamental pair i, j with Tj“, leading to the well-known 
formula II43[I44II 


jyjar-an ^ ^ Tr{T^-(')T“-m • • • T“-(-i)T^-} A<,. (15) 

(JGSn-l 

The formulas above are not quite the same as in the mentioned references. One ingredient that 
is missing is that, if all gluons are on-shell, the partial amplitudes are given by a single 
amplitude with permuted arguments: 

Aff(l,...,n — l,n) = A(ct(1 ),... , cr(rL — 1 ),n) . (16) 

Each argument in the form of a number, say I, refers to both the momentum and the helicity of 
gluon 1. This is easy to understand considering that all external gluons are essentially equivalent. 
If some of them are off-shell, this does not seem so obvious, but we may consider the case in 
which all of them are off-shell, and thus equivalent. Now, each number I refers to the momentum 
and direction, or longitudinal momentum component, associated with off-shell gluon 1. We may 
take the on-shell limit of gluon 1 , which in ifSTll is argued to lead to 

^y^ar-an ^ Tr {T“-n) 7^2) ... yi(a(icT(n - 1 ),n) 

O-eSn-l 

+ Y_ Tr{T“-n)T^a(2) ...Ta.(n-i)Tan| (17) 

O-eSn-l 

It is explained in ifSTll how this coherent sum over helicities turns into an incoherent sum for the 
squared amplitude. We see that for each helicity, the color decompostion including Eq. (fT^ . still 
holds. The color decompostion presented in [[39ll also holds, since it is proven there using only 
group theory arguments, and is independent of the exact form of the partial amplitudes. 
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4 Usage 

The program can be obtained from 

http://bitbucket.org/hameren/amp4hef 

The main directory contains a readme file describing how to compile and use the program. There 
is an example directory with a simple Monte Carlo program to show how to use the program 
explicitly. The program is written in Fortran 2003, but it eventually provides a module that 
makes a few subroutines available that only take (arrays of) intrinsic integer, real and complex 
type variables as arguments. More specifically, the 

module amp4hef 

provides first all the 

subroutine init_amp4hef 

which does not take any arguments, and must be called once before any other routine. Further¬ 
more, it provides the 

integer,parameter :: fltknd=kind{IdO) 

which is set to this value in the module amp4hef_spinors, which can be found in the source 
file amp4hef .spinors . f 90. All real and complex variables are of this kind. Then, the module 
provides 

subroutine construct_ng{ Noffshel ,Nonshell ,momenta ,directions ) 
integer ,intent (in) :: Noffshell ,Nonshell 

real(fltknd),intent(in) :: momenta{0:3,*) ,directions{0:3,*) 

This routine takes external momenta and directions as arguments, and prepares all spinors etc. 
for the evaluation of amplitudes and matrix elements. The input Noffshell refers to the 
number of off-shell momenta, and should be 0, 1 or 2. The input Nonshell refers to the 
number of on-shell momenta. All momenta should be provided via the array momenta. The 
size of the second dimension should be at least n=Nof f shelH-Nonshell. The initial-state 
momenta should be momenta {p., i) =kf for ie{1,2}, while the final-state momenta should be 
momenta {p., 1) =pf for iG{3,..., n}, where p.G{0,1,2,3}. The size of the second dimension of 
the input directions should be at least Noffshell. This array should provide the directions 
associated with the off-shell momenta: directions {p., 1) They should satisfy Eq. ([8]). 

After calling subroutine construct_ng, the squared amplitude summed over colors and 
helicities can be evaluated for the given momenta and directions with 

subroutine matrix_element_ng{ rslt ) 
real{fltknd),intent(out) :: rslt 

The output of this routine is missing a factor where n is the total number of gluons. 

The number of colors is fixed to Nc = 3. This can be changed in the file amp4hef_ng. f90. 
The output of this routine is normalized such that when kf ^ pf for each off-shell gluon, then 
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the result beeomes the standard result with on-shell gluons only, including the sum over all their 
helicities. 

Also after calling subroutine construct_ng, color-ordered helicity amplitudes can be 
evaluated with 

subroutine amplitude_ng{ rslt ,helicity ,permutation ) 
complex(fltknd),intent(out) :: rslt 
integer ,intent{in) :: helicity{*) 

integer ,intent{in) :: permutation{*) 

The size of the array helicities should be at least n^Nof f shell+Nonshell. Its entries 
refer to the helicities of the on-shell gluons and follow the same enumeration as the momenta 
given to subroutine construct_ng. Their values should be one of—1 or-|-l. The size of the 
array permutation should be at least n — 2, and only refers to the final-state gluons. Denoting 
Ai=helicities (1) for iG{1 ,2,..., n}, and cr(i)=permutation (i) for iG{1 ,2,..., n — 2}, 
this routine returns the value of 


P2+o-(l) , ^2+o-(l) ; P2+o-(2) ) A2 +o-( 2) '> • • • i P2+o-(n-2) ) A2+o-(n-2) ) • (18) 

Besides the square root of the factor mentioned before, the output of this routine is also missing 
a factor \/\kl\ for each off-shell gluon. 

Similarly to the routines described above, the module amp4hef also provides 

subroutine construct_ng_xpr(Noffshel,Nonshell,momenta,directions) 
subroutine matrix_element_ng_xpr(rslt) 

which are defined in the file amp4hef_ng_xpr. f 90, and which use analytic expressions rather 
than numerical BCFW recursion. The total number of gluons is restricted to a maximum of 5 for 
these. To give an impression of the size of the expressions, we present here one of them: 


Air, 2% 3-, 4-, 5+) 

1 


1 [12] ^(43) 3 

" Ki K2 (5|Ki +K 2 I 2 ](3|Ki +^211](54) +]t 2 V 

( 32 )^[ 51]3 1 [ 25 ]'^( 21)3 


KiK| (2|K2+?S3|4](3|K2|1](K2+f^3)2 (1 |Ki+K 2 I 2 ](2|Ki|5](2|Ki|2][23][34][45] 

J__ 

K2 (l|Ki+K2|2][2Kizi3+lzi4)(Ki+?S5)|1](5||i3+f^4|2](15)[23][34]((1|K2|1][24] + (13)[34][21]) 


1 [51]3(2|^4+?S3|2]3 

+ Ki (2|K2+f<3|4](2|Ki|5][2K^4+T^3)(Ki+?S5]|1][23][34](Ki+f<5)2 

1_[12P(13)-^[51]^_ 

K2 (3|K2n](1IK2|l](3|Ki+K2n](llK2+?S3n][45]((l|K2|l][24] + (13)[34][21]) 


When both off-shell gluons go on-shell, the first term corresponds to the helicities (1+,2‘''), 
the second to (!■*■, 2~], while the third and fourth combine to [1“,2^). The last two terms are 
irrrelevant in that limit, and the helicities (1“,2“) correspond to a vanishing amplitude in the 
on-shell case. 
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The module amp4hef also provides the routines defined in the file amp4hef_ons_xpr. f90, 
whieh eolleets some expressions from literature for squared amplitudes summed over eolors and 
helieities [|45ll . These are mainly ineluded to elarify the normalization eonventions used in this 
program. 


5 Validation 

The Monte Carlo program in the example direetory simulates the generation of events by read¬ 
ing event files that were produeed with the help of A Very Handy LIBrarjll] in Fortran. The files 
inelude the event weights obtained with AVHLIB, and have been used to validate the squared am¬ 
plitudes summed over spins and eolor provided by AMP4HEF. AVHLIB provides a eompletely 
independent implementation for the numerieal evaluation of seattering amplitudes, employing 
Dyson-Sehwinger reeursion, for arbitrary proeesses, both with on-shell and off-shell initial-state 
partons, and ineluding both QCD and eleetro-weak interaetions. Furthermore, it provides an 
effieient phase spaee generator II46[I47H and other tools for the evaluation of distributions of arbi¬ 
trary observables using the Monte Carlo method of integration, with whieh the aforementioned 
event files were produeed. The amplitudes ean be provided in several representations by AVH¬ 
LIB, among whieh is the deeomposition into eolor-ordered helieity amplitudes. So besides the 
squared amplitudes, also the eolor-ordered helieity amplitudes from AMP4HEL were validated 
seperately with the help of AVHLIB. 

6 Summary 

We presented a program to numerieally evaluate multi-gluon amplitudes, whieh ean be applied 
in faetorized eross seetion ealeulations for hadron eollisions whieh ask for off-shell initial-state 
gluons. It ealeulates squared amplitudes summed over eolors and helieities for up to six exter¬ 
nal gluons with up to two of them off-shell. Analytieal expressions for eolor ordered helieity 
amplitudes are employed for up to five external gluons with up to two of them off-shell. Higher 
multiplieities are evaluated using numerieal BCLW reeursion. The program has been validated 
with an independent library employing numerieal Dyson-Sehwinger reeursion. 
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A Construction of the direction 


Given the Sudakov deeomposition in terms of two light-like momenta j >2 

k*^ = Xi p!|*^ + X 2 P 2 + kj with pi-ky = p 2 -kT = 0 and pi ■p 2 = ^ s 7 ^ 0 , ( 20 ) 

we may ehoose 

P‘‘=pr-^p;-xK wi,h . ( 21 ) 

We then have 

= xip^^ + k;^ ( 22 ) 

with p^ defined above, and with 

k|^ = (1 + xxi ] k^;^ + ^ ^<^1 + ^<^ 2 ^ P 2 ) (23) 

satisfying 

p-p=p-k|=-0. (24) 

Notiee that x —> 0 for X 2 —> 0, so the eonstruetion is eonsistent. Also, if k^ beeomes light-like, 
that is if X 1 X 2 S -|- kj —> 0 , we have x —> —1 /xi and k|^ —> 0 . 

Alternatively, if k^ is given as the sum of two light-like momenta. 




(25) 


we may ehoose 

p^ = p^ - P 2 - ^ (pi ly^lpJ + j (p2ly^lpi] (26) 

for any z 7 ^ 0 , or in terms of spinors: 

Ip) = Ipi)+z|p2) , Ip]-= Ipil - - Ipi] . (27) 

z 

A deeomposition like Eq. (|25]) ean always be eonstrueted given any auxiliary light-like momen¬ 
tum with q-k 7 ^ 0 : 
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